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We calculate the imaginary part of the retarded two-loop self-energy of a static vector boson 
in a plasma of quarks and gluons of temperature T, using the imaginary time formalism. We 
recombine the various cuts of the self-energy to generate physical processes. We demonstrate how 
cuts containing loops may be reinterpreted in terms of interference between 0(q) tree diagrams and 
the Born term along with spectators from the medium. We apply our results to the rate of dilepton 
- - - , production in the limit of dilepton invariant mass E » T. We find that all infrared and collinear 

■ singularities cancel in the final result obtained in this limit. 
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. I. INTRODUCTION 

oo 

(N : The imaginary parts of retarded self-energies represent extremely important quantities in thermal field theory. They 
provide information about various quantities of physical interest in the medium. Primary among these are the decay 
^ • and formation rates of particles Boson self-energies provide information about quantities like Z decay rates 

and production rates of dileptons and real photons Q from a quark-giuon-plasma (QGP). The spectrum of lepton 

■ pairs {i.e., e+e^, and real photons emanating from such a plasma has been considered as a promising signature 
of QGP formation [QH]. This owes to the fact that the photons or dileptons suffer essentially no final state interaction. 

Some years ago the contribution to the rate of dileptons produced at rest in the plasma at first order in the strong 
coupling constant, was evaluated ||^. This included reactions like three particle fusion {qqg — * 7*), Compton scattering 
(qg — > qj* or qg — > ^7*), pair annihilation (qq —>■ gj*), Born term with vertex correction, and Born term with quark 
or antiquark self-energy correction. This calculation was performed in the real time formalism, both in a Feynman 
' diagram approach in thermo- field dynamics, and by taking the imaginary part of the two loop photon self-energy. In 
^pH, the case of massless QCD, each of the contributions mentioned above contain infrared or collinear singularities. These 
were regulated at intermediate stages of the calculation by giving masses to the quarks and gluons. The combined 
D , rate from all these processes was then found to be free of all divergences in the limit of vanishing masses. This 
1-^ calculation was also performed simultaneously by another group Q, who dimensionally regularized the singularities 
J> [ at intermediate stages of the calculation. The end result remained the same: when all the different processes were 
• '~j ' summed, the divergences cancelled and dilepton rate at next-to-leading order remained finite. 

rxj Recent calculations employing a multiple scattering expansion, however, have found a remnant collinear divergence 

[^,||. This result has been commented upon and the issue of divergences remained unresolved In the wake 
■ ■ ■ ' of this strife, we revisit this problem in a systematic calculation. Also, to the best of our knowledge, a complete 
calculation of the imaginary part of a heavy vector boson retarded self-energy in the imaginary time formalism has 
yet to be performed. This is the subject of this paper. The scalar boson self-energy was examined recently ]lT[ . There 
are various advantages to such a calculation: the basic Feynman rules are easily generalized from zero-temperature; 
there is no doubling of degrees of freedom and no matrix structure of propagators; multiple poles which lead to ill- 
defined products of delta functions in the real-time formalism are easily and naturally handled both in the Matsubara 
sums and in the analytic continuation. The purpose of this calculation is thus many-fold. A first goal is to enumerate 
and interpret the various physical contributions contained in the imaginary part of the two-loop self-energies. In doing 
this, it shall then be shown that cuts containing loops may be re-expressed as interference between tree diagrams and 
Born term with a thermal medium spectator. Importantly, we also demonstrate how double poles may be simply and 
elegantly dealt with, in the Matsubara sum and in the analytic continuation to real energies. We finally concentrate 
on eventual collinear and infrared divergences in the ensuing rates. In this study, we focus on the singularity structure 
in the region of phase space investigated by the authors of Refs. |^,^. Even though we explicitly calculate the 
self-energies of static virtual photons, the results may be easily applied to other vector bosons in-medium, with the 
exception of the gluon which admits other self-energies in a QGP. 

The various sections are organized as follows: in section II we begin by evaluating one of the self-energy diagrams 
of a static photon with an imaginary energy at two loops (the impatient reader may skip ahead to section VI where 
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the various cuts of the self-energy are recombined to provide physical interpretations of the various terms obtained; 
following which the infra-red behaviour of heavy photon production will be discussed); in section III we analytically 
continue the photon energy to real values and obtain the imaginary part of the corresponding retarded self-energy; 
in section IV we evaluate the other self-energy topology; in section V we analytically continue this self-energy to real 
values of photon energy and find the retarded imaginary self-energy; in section VI we combine the tree-like cuts and 
reinterpret them as physical processes with thermal distributions on the phase space factors; in section VII we attempt 
to interpret cuts containing loops in terms of the recently proposed spectator interpretation Jl2| ; in section VIII we 
take the limit of heavy-photon production (E >> T), and evaluate the various contributions; in section IX we combine 
all cuts, demonstrate the cancellation of the collinear and infrared divergences, and present our results; we present our 
conclusions and brief discussions in section X. A few short appendices follow. In the interest of quantitative accuracy 
and repeatability, we have presented many calculational details: the issue being addressed here is technical, and thus 
demands a rigorous treatment. 



II. THE SELF-ENERGY: TOPOLOGY I 




FIG. 1. The first topology for the self-energy. 



We evaluate the photon self energy with a gluon running across as shown in Fig. |l|. To begin with, we derive the 
expression for the effective quark photon vertex corrected by a gluon running across, i.e.. 

Where e may be taken to be the electric charge of the quark. In standard notation, the expression for the effective 
vertex in Feynman gauge may be written down as 



pi. J (2^)3 ,2 i^k^ain _ ^ _ i^e^nj^i^h,,91 ). (2.1) 



The Matsubara sum in the effective vertex may be simply evaluated using the method of Pisarski |13|. In our 
notation (see Appendix A, see also ||l^), this is given in the static limit (p = 0) as, 

4 J (2^)3 ^^^^^^ qE,^,E,^k{p" - is2 - S3)E,^k) 

j {si + S2)/2 - sinjEg^k) + S2n{q) (si + S3)/2 - sin{Eg_k) + Ssnjq) \ 

kO - s,q - S2Eg.k kO-pO-s,q-ssEg., /• ^'-'^ 

Where Si, 82,83 are sign factors which are summed over the values of ±1. We may now use the above result to 
write the full self-energy of the photon in the static limit as 
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/'tfJ (2') 



(2.3) 



Where, kg stands for the four component quantity: 



We note that the effective vertex may be written as 

to highhght the structure of 7 matrices contained within it. The trace of the 7 matrices is straightforward. This 
gives the fuU self energy as 



feO 



(fik 



H I (27r)3 ^ 



(2.4) 



For convenience we change S2 — > —S3 and S3 — > — S2- To evaluate the Matsubara sum we follow the method of 
reference lid. This method converts the Matsubara sum into a contour integration in the complex plain of fc*^, i.e., 



2_2 fioo+n 



2ni 



-[l/2-n(fc-)]- 



S5 



(g-fc),,(g-fc),3 

g(P° - (S2 - S3)Eq-k) 



S2- 



(si - S3)/2 - sin{Eq^k) - S3n{q) 



ssfcO - siq + s^Eq-k 



+ S3 



(si - S2)/2 - sin(Eq^k) - S2n{q) 



p" - S5fc" + Sig - S2Eq^k 

Where, {q — k)^ stands for the four component quantity: 



S4fc/3,; 



ssfcO — — Sik 



(2.5) 



qx kx qy ky qz k^ 



q sm t) cos < 



q sm sm ( 



q cos 6 — k 



\q-k\' \q-k\' \q - k\ 



^/W+(f^2kqcose ^A:2 + g2 _ 2kqcoa8 y^k'^ + q^ - ikqcosO 



The A:*^ integration is from ~oo — > 00 on the positive side of the real axis. We may thus close the contour on the 
positive side. Note that the function is vanishing as — > 00. The result of this integration will simply be the sum of 
the residues at the corresponding poles. Looking at the above expression we note that the pole structure is different 
depending on whether the term being considered is the first or the second one in the curly brackets. We note the 
following poles: 

i) 1st order pole: at fe" = fc, requires S5 = s (in both terms). 

ii) 1st order pole: at A;" = A; + ssp", requires S5 = S4 (in both terms). 

iii) 1st order pole: at fc" = sssig — s^ssEq-k, requires sssig — s^s^Eq-k > (only in the first term). 

iv) 1st order pole: at fc" = ssp*^ — szS2Eq^k + sssiq (only in the second term). 

In the following, each of the poles are evaluated in a separate subsection and then summed up. 
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A. 1st order pole at fc" = fc 



This is the pole of the first outer propagator {i.e., not a propagator in the effective vertex), it is a pole for the entire 
self energy expression, ft has the obvious residue of 

d^kd^q, 



n{:(^)=4eV/^[l/2-n(fc)]f^ 



(g-fc).,(g-fc),3 

q{P° - (S2 - S3)Eq^k) 



32- 



(si - S3)/2 - sih{Eg^k) - S3n{q) 



sk - siq + s^Eq^k 



S3 



jsi - S2)/2 - sin(Eq-k) - S2n{q) 
- sk + siq - S2Eq^k 



.sk — — s^k 



(2.6) 



Note that there is an extra negative sign in the residue as the contour is being taken in the clockwise sense. 



B. 1st order pole at k° = k + s^p^ 

This is the pole of the second outer propagator, it is a pole for the entire self energy expression. It gives the residue, 

d^kd^q ,^ _ ska s 



n;i(i?) = 4eV/^[l/2-n(fc)]-, 



32 



Sik — sk 
(si - S3)/2 - sin{Eq^k) - szn(q) 



S4,k - siq + SsEq-k 



S3- 



(si - S2)/2 - sin{Eq-k) - S2n{q) 



-Sik + siq - S2E1 



q—k 



s^k 



S4 



(2.7) 



Note that the in the distribution function has been dropped. This may be done as e^"^ = 1 {p^ is a discrete 
even frequency), and, secondly, as we are eventually going to analytically continue the self-energy to complex values 
of The correct analytic continuation is given by that function which has no non- analytic behaviour off the real 
axis ]l6[ |. One may easily check that the above function with a in the distribution function will have poles at 
p° = —k + i2{n + 1)ttT. Note that in this pole we may switch 

S — > -S4, S4 -S, S2 -S3, S3 -S2, Si -Si, 

and noting that k^s{q — = ks{q — k)^^, we find 



C. 1st order pole at k" = s^siq — S5S3E, 



q — k 



In the expression for the effective vertex (Eq. (2.2) or Eq. (2.5)), we note the presence of two terms inside the curly 
brackets with different pole structures. This is a pole of the first term, and is realized only if s^siq — s^s^Eq-k > 0. 
Thus the residue is 



-^^-^-[1/2 - n{s5Siq - ssSsEq-k)] 



Skn 



siq - ssEq^k - sk 



(g-fc)s.(g-fc).3 

g(p° - (S2 - S3)Eq^k) ' 



S2- 



(si - S3)/2 - sin{Eq_k) - s^niq) 



S5 



S4fc/3,S4Q(s5Slg - S5S3£'g-fc) 

sig - .s^Eq-k -p" - S4fc 



(2.8) 
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D. 1st order pole at fc" — ssp" + sr,siq — S5S2-B, 



q — k 



This is the pole of the second term in the curly bracket mentioned in the previous section. It is realized only if 
SbSiq — SQS2Eq^k > 0. The residue is 



/d'^kd'^a 



q - S2Eq^k 



sk 



iQ-k),(q-k), 



S3 



(si - S2)/2 - Sin{Eq_k) - S2n{q) 



qiP° - (S2 - S3)Eq^k) 

Note that in this pole we may switch 
S2 -S3, S3 
With this operation, we find 



-S2, Si 



-S5 



-si, S5 



Sikfj^siQjssSiq - S5S2Eq^k) 
siq - S2Eq^k - Sik 



-S5, s 



-S4, S4 



(2.9) 



Thus the full photon self-energy to second order in the coupling constant for the diagram of Fig. [l] is given by 
summing up the results of the preceding four subsections, i.e., 



III. IMAGINARY PART OF THE FIRST SELF-ENERGY TOPOLOGY 

We now proceed with evaluating the discontinuity in the first self-energy as is analytically continued towards 
the positive real axis from above i.e., E + ie. Analytically continuing p'^ will give us the retarded self energy of 

the photon in real time in terms of a real continuous energy p° = E. The expressions to be continued are nj^(^) and 
n^(I?). The presemce of the theta function in n|^(-D) complicates the pole sturcture that one would obtain during 
analytic continuation of this expression. Using standard techniques (outlined in Sect. IV. D.), we decompose the 
theta function to obtain. 



n;i(z)) = -4eV/ 



d^kd^q SS4S3 



(27r)6 q 

\ [ks ■ {q^k)^J[k,, ■ (g'^fc),J[l - njEq^k) + n{q)][l/2 - h{Eq_k + g)] 
I [pO - {sk - s^Eq-k ~ S5q)][p0 - (-S5 - S3)Eq-k\[s5{Eq-k + q) - Sik] 

[ks ■ {q^k)J[k, ■ {q'^k)^J[h{Eq^k) + n(g)][l/2 - n(g - g,,,.)] \ 
[pO - (sk + ssEq^k - S5q)][p° - (55 - S3)-Bq_fe] [35(9 - Eq^k) - S4A;] J " 

Analyzing the expressions for 11^ (A) and nj;J(£'), we note the following discontinuities: 
Poles of type p° = 2k: 

i) 1st order pole in 11^ (A): at p'~' = 2k, requires s — — S4 = 1 (in both terms that make up n|^(A)). 
Poles of type p° = 2Eq-k- 

ii) 1st order pole in 11^ (A): at p^ = 2Eq^k, requires S2 = —S3 = 1 (in both terms that make up 11^ (A)). 

iii) 1st order pole in n^(_D): at p° = 2Eq^k, requires S5 = S3 = — 1 (only in first term of 11^ (£>)). 

iv) 1st order pole in ^{D): at = 2£'q_fc, requires S5 = —S3 = 1 (only in second term of ^{D)). 



(3.1) 
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Poles of type p" = sk + siq + S2Ek-q: 



v) 1st order pole in n|^(A): at p° = sk — Siq + S2-Bg_fe, requires sk — Siq + S2-Bg-fe > (only in the second term of 

vi) 1st order pole in 11^ at — sk — s^q — Sc,Eq_f^, requires sk — s^q — s^Eq_k > (only in the first term of 

n$l(s)). 

vii) 1st order pole in n^(i?): at p'^ = sk — s^q + s^Eg-k-, requires sk — s^q + s^Eg-k > (only in the second term 

ofn^i(B)). 

We may write down the expression for 2Il'^{A) highlighting its real and imaginary parts as p° — > -E + ie as 

d^kd^q S4[ks ■ {q^k)J[k, ■ {q^k)J 



(27r)6 
1 



X < S2 



E - {s2- S3)Eg_k 

{si - S3)/2 - sifijEq^k) - S3n{q) 
sk - siq + ssEq-k 

1 

E - sk + siq - S2Eq^k 
1 



Eq+k 

iTr5{E - (s2 - S3)Eg_k) 



+ S3((si - S2)/2 - sin{Eg_k) - S2n{q)) 



- m5{E — sk + siq — S2Eq-k) 



>[l/2 -n{k)] 



E — {s — Si)k 



— i'ir6{E — {s — S4:)k) 



(3.2) 



We now write down the various discontinuities as enumerated above (Note: we are now looking at 211^ (A) and 
21I'j^{D), so the overall factors have doubled): 



Disc[2U';^{A)]a = {+27Ti)8e^g^ J 



cfkdS ~2[k+-{q-k)J[k.-{q-k)^ 

{2Tlf q(E - (S2 - S3)Eg-k) 



^^^ (5i-.3)/2-.in(i^, 0-.3n(g) _ _ 
k — siq + ssEq-k 



n- roTTMMM (^0 2 2 f d^kd^qS4[ks-{q-k)^][h,-{q-k)_] 
Disc[2U'i^{A)\b = (+27r«)8e g J ^ 

-SiE-2Ek-g). (3.4) 



(2^) 

[I ~ 2fiiEg_k)]2q I [l/2-n(fc)] 
{sk - Eq_k)^ - q^ \E-{s-S4)k' 



n- [oTTMmM r^o 2 2 ffkfql [fc. ■ (g - fc)^][fc, ■ (g - fc)_] 
^'''^^Km^ = (+2-)8e , J - I El_,-{sk + qY 

[ks-{q-k)^][ks-{q-k)^]\( / mm /o -it? , m 

+ q--{Eq.k-skY / 1 - + "(^^^ f'/' - "(^^-'^ + 

- [n(£;,_fc) + n(9)][l/2 - n{q - Eq_k)]}5{E - 2Eg_k)- (3.5) 

We now proceed to the evaluations of the discontinuities of the type = sfc + S2-Eg-fe + siq. We first change 
Si — > — Si in 2nj|^(A), and S5 — > — S5 in the first part of 2nj|^(D). Hence the discontinuity occurs in 2nj|^(A) at 
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p° = sfc + S2Eq-k + siQ only when sk + S2Eq^k + siq > 0. This may happen in only one of four instances: 



the discontinuity occurs in the first term when s = +,S5 = + or when s = — ,55 = +; in the second term when 
s — S5 — + or when s = +, S5 = — . 

Thus the discontinuity in 2nj^(A) occurs in four parts: 



Disc[2ni/^{A)]c = (+27rj)8eVy" 



^[l/2-n(fc)] 



{2TTfq 

S3S4[fc+ ■ (g - k)^][ks^ ■ [q - fc)^J[-l + n(£'q_fc) - 

S3S4[fc- • [q - k)^][ks^ ■ {q - k)^J[-l + n{Eq-k) - n{q)\ 
[-k + q + s^Eq-k] [q + Eq^k + s^k] 

s3S4[fc+ ■ (g - fc)_][fcs4 • (g - fc)s3]["-(-^g-fc) + 

[fc + g + s^Eq^k\[-Eq^k +q + Sik] 
S3S4[fc+ ■ [q - fc)+]fe4 ■ (g - k)^^][-n{Eq.k) - »(g)] 

[fc - g + SzEq^k\[Eq-k -q + S4:k] 



5{E-k- Eq-k - q) 

6{E + k- Eq^k ~ g) 

6{E^k + Eq^k - q) 
5{E-k-Eq_k+q)\. 



The discontinuity in 2n'J(I?) also occurs in four parts: 



Disc[2m{D)]c = (-271^)86^5^ 



(27r)6g 



[1/2 - n{Eq-k + q)] 



\ szSi[k+ ■ {q - k)^][ks^ ■ (g - fc)^J[l - njEq^k) + n{q)] 
I [k + q + ssEq^k] [g + Eq^k + S4k] 

_ S3S4[k- ■ (g - fc)^][fcs4 ■ (g - - n{Eq^k) +n{q)] 

[-k + q + S'iEq^kWq + Eq^k + Sik] 
+ [1/2 -h{q-Eq^k)] 

j S3S4[fc+ • (g - k)_][ks^ ■ {q-k)^J[h{Eq_k)+n{q)] 



6{E-k- Eq^k - g) 
5{E + k- Eq_k - q) 



[fc + g + S:iEq-k][-Eq-k + g + s^k] 

S3S4[fc+ ■ (g - fc)+][fcs4 ■ (g - k) ^^\[n{Eq_k) + n{q)\ 
[k-q + S'iEq^k] [Eq^k - g + Sik] 



6{E-k + Eq_k - g) 



6iE ~k~ Eq_k + g) 



(3.6) 



(3.7) 



IV. THE SELF-ENERGY: TOPOLOGY II 

We begin by evaluating the photon self energy with one quark line containing a gluon loop. The quark self energy 
may be written as: 

- ,m,. - E / (4-1) 

Using the identity 7^ ~ —2 q being the momentum of the gluon, the Matsubara sum in the quark self energy 
is evaluated using the method of Pisarski [|l3| . It is given in our notation( Appendix A) as 

9%,tlk5'''' f A {^fi)s.. jsi + S2)/2 - s,n,{Ek-q) + ^2^(9) . 

^ 2 i (^ ^ —q P - s^Ek-q - s,q ^ ^^-2) 

S1S2 ^ 

Where, the self energy has been written in the final form to highlight its matrix structure. We may use this to 
write the full self energy of the photon in the static limit as 
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FIG. 2. The second topology for the self-energy. 



13^ J (2^)3^ ' ^ 

' k° ^ ' S1S2S3S4S5 L 



67^ 



2(fcO - Safe) 



67'^ 



2(/fcO-S4A:) 2(/fcO-pO-S5fc) 



(4.3) 



We choose the z direction to be defined by the direction of k. Note that kg^ ■ ks^ = —2Ss^,-s5- Note also that the 54 
and S3 dependence of the photon self-energy is identical. This allows us to write down the self-energy as, 



feO 



_fk_ 

(27r)3 



2e' 



2(-S3)Sjj ■ ks 



2( S5)'^jj ' 



(fcO - S3fc)(fc0 - Sik){kO - pO + S4fc) (fco -pO - S5A:)((fcO)2 - fc2) 



(4.4) 



Where summation is implied over all the sign variables present i.e., si, S2, S3, S4, S5. Note that the double pole is only 
present in the first term. 

We now need to evaluate the Matsubara sum over fc°. For this we follow the method of reference This method 
converts the Matsubara sum into a contour integration in the complex plain of fc". The color factor from the quark 
self energy combined with that from the rest of the diagram becomes tr[i'^, t^]5°'^ — 4. Using this we obtain the self 
energy of the photon as: 



^eV n"-' f d'kd^l [1/2 - n(fc°)][(5i + S2)/2 - s^Ek^g) + •^2^(9)] 



27ri 



q[sk° - S2Ek^q ~ siq] 

Sb[k-q)^^ - ks, 

[sk° - S3k][sk° - S4)fc][sfc0 -pO + S4A:] [sk"-p" - S5fc][(fc0)2 



-ioo+e J (27r)6 



(4.5) 



The k^ integration is from — cx) ^ 00 on the positive side of the real axis. We may thus close the contour on the 
positive side. Note that the function is vanishing as 00. The result of this integration will simply be the sum of 

the residues at the corresponding poles. Looking at the above expression we note the following poles: 

i) 2nd order pole: at k^ — k, requires S3 = S4 = s (only in the first term). 



ii) 1st order pole: at k^ — k, no requirement (only in the second term); requires S3 — — S4 
(only in the first term). 



s or —S3 



Si 



iii) 1st order pole: at fc" = fc -F sp", requires S4 = — s (only in first term); requires S5 = s (only in second term) 

iv) 1st order pole: at k^ = ss2Ek-q + ssiq, requires ss2Ek-q + ssiq > 0. 

In the following each of these poles will be evaluated in a separate subsection and then summed up. 
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A. 2nd order pole at k° = k 



We begin by evaluating the 2nd order pole. The origin of this pole can be traced back to the two propagators 

whick may go on-shcll sinmltancously In the real time formalism this leads to the ill-defined sqnare of the Dirac delta 
function. In imaginary time, however, this pole is easily dealt with: the residue of a function f{k^) at a second order 
pole at fc'^ = A; is simply given as -;^{k^ — k)^f{k^)\k°=k- Using this we get the residue of IT at this pole as 



TJt^(A\ _ S,2„2 f d^k(fl [(51 + S2)/2 - Sih{Ek-q) + S2n{q)] 



s{l/2-n{k)y 



l/2-n(fc) 



[sk - S2Ek-q - siq] [pO - 2sk] [sk - S2Ek-q - siq]'^\p° - 2sk] 
l/2-n(fc) 1 



[sk - S2Ek-g - siq]\p° - 2sA;]2 



(4.6) 



Note that as in the case of the first self-energy topology there is an extra negative sign in the residue as the countour 
is closed in the clockwise sense. 

B. 1st order pole at fe" = fe 

This obvious residue may be easily evaluated using the methods outlined in the previous section. 



WJB) = -8eV 



d^kd^q [(si + S2)/2 - sin{Ek-q) + s2n{q)\[l/2 - n{k)] 



(27r)6 



q[sk - S2Ek-q - siq] 



-sjk - q)^^ ■ ks s{k - q)^^ ■ k-s S5{k - q)^^ ■ ks. 



1 2fcp0 2fc[pO-2sA;] 2k[p" - (s - S5)k] j' 

In the above, we sum over the two possibilities of S5 = ±s to get the factor in the bracket as 



(4.7) 



-s{k - q)^^ ■ ks s{k - q)^^ ■ k-s s{k - q)^^ ■ ks 

T ^ 7 r n ^ J 1 "r 



2kpO 



2k\p° - 2sk] 2kp° 
-s{k - q)„^ ■ k^s \ 



+- 



2fc[p0 - 2sk] 



Hence, m{B) = 0. 



C. 1st order pole at fc" = A: + sp° 



This gives the residue 



n^rr^ = -R^ ^ f ^'fc^'g K^i + '^2)72 - sMEk-g) + S2n{q)][l/2 - n{k)] 
^[^) ^eg J q[pO-{s2Ek-q + Siq-sk)]\pO + 2sk] 

j S3(fc-g),, - ks ik-q)s2 -k] 



1 s{p° - (s3 - s)k) p' 
Switching s — > —s, and summing over S3 = ±s we get. 



2 2 / d'^kd^q[{si+S2)/2-sin{Ek-q) + S2n{q)][l/2-h{k)] 



(27r)6 

{k - q)s^ ■ k 
{pO - 2sk) 



- {s2Eu-q + siq + sk)]\p^ - 2sk] 



(4.8) 



(4.9) 
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D. 1st order pole at fc" = ss2Ek^q + ss\q 



This pole is realized only if ss2Ek-q + ssiq > 0. This condition may be enforced with the following set of delta and 
theta functions: 

Ss,S2Ss,s, + Ss,s2Ss.-s,e){Ek-g -q)+ Ss,-s2Ss,s,ei{q - Ek-q). (4.10) 

We start with the second and third terms: 

(Pkd^q [1/2 - - q)][n{Eu-q) + n{q)] 



(27r)6 q 



[sEk-q - sq - s:ik][sEk-q - sq - s^k][p° - {s^k + sEk-q - sq)] 
S5{k - q)^ ■ ks^ "I 



[pO - {sEk-q - sq - s^kMEk-q - qf - fc2 ^ 
Similarly we find for the third term 

d^kdS [1/2 - fi{q - Ek-q)][~n{Ek-q) - n{q)] 



Q{Eu-q-q). (4.11) 



n;i(i^,3) = 8eVy ■ 



(27r)6 q 



[-sEk-q + sq - s:ik][~sEk-q + sq - S4k][p" - {sik - sEk-q + sq)] 
S5{k - g)_3 • fcss 



[pO - i-sEk-q + sq- s^kMEk-q - qY - A:2 



Q{q-Ek-q). (4.12) 



Now, switch s — > — s in the 3rd term, and noting that 1/2 — n{q — Ek-q) = — 1/2 + n{Ek-q — q), we note that the 
second and third terms can be combined to give 



n;i(i?,2 + 3) = 8eVy 



d3fcd3g [1/2 - hjEk^q - q)\\h(Eu-q) + n(g)] 
(27r)6 q 

sz{k - g)^ ■ fc^3 

\sEk-q ~ sq ~ S3k][sEk-q - sq - Sik]\p° - (s^k + sEk-q - sq)] 

55 (fc - q)s ■ ks, 1 

' [pO - {sEk-q - sq - S5k)][{Ek-q - g)2 - fc2] f • 



(4.13) 



Note the absence of the theta fu nction s in the above equation. Now we may also write down the residue from the 
first set of delta functions in Eq. (4.10) as. 



U^^iD,l) = 8e^9' I 



<fikd\ [1/2 - n{Ek-q + q)][l - njEk-g) + n{q)] 
(27r)6 q 

ssjk - q), ■ fcsa 

[sEk-q + sq - s^k][sEk-q + sq - S4,k][p" - {s^k + sEk-q + sq)] 

s^jk - q)^ ■ ks^^ 1 

' [pO - {sEk-q + sq- s,k)][{Ek-q + g)2 - f ■ 



(4.14) 



The total expression obtained by summing up the results from the preceding 4 subsections will give us the full 
self-energy of the photon to second order in the coupling constant for the diagram of Fig. ||, i.e., 

n;i = w^iA) + u^JB) + u^JC) + u^;iD, i) + ^d, 2 + 3). 
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V. IMAGINARY PART OF THE SECOND SELF-ENERGY TOPOLOGY 



We now proceed with evaluating the discontinuity in the second self-energy as p° is analytically continued to a 
positive real value i.e., — > E + ie. Analyzing the expressions derived in the above sections we note the following 
discontinuities: 



a) Poles of type p° = 2k: 

i) 1st order pole in nj[^(74): at p° = 2k, requires s = 1 (only in first and second terms). 

ii) 2st order pole in n^(A): at p° = 2k. This occurs in the third term in the bracket and requires s = 1. 

iii) 2nd order pole in n^(C): at = 2k, requires s = 1 and S3 = 1 (only in first term). 

b) Poles of type = sk + siq + S2Ek-g: 

iv) 1st order pole in nj^(C): at p^ = sk + siq + s^Ek-q, requires s = Si = S2 = 1, or — s = si = S2 = 1, or 
s = —Si = S2 = 1, or s = Si = — S2 = 1 (in both terms ). 

v) 1st order pole in njJ(-D, 2 + 3): at p° = 54^ + sEk-g — sq, requires S4 = s = 1 or S4 = — s = 1 (only in first term); 
at p^ = —s^k + sEk-q — sq, requires — S5 = s = lors5 = s = — 1 (only in second term). 

vi) 1st order pole in nj|^(-D, 1): at j*" = S4k + sE^-q + sq, requires S4 = s = 1 or — S4 = s = 1 (only in first term), 
at p*^ = —s^k + sEk-q + sq, requires — S5 = s = lors5 = s = l (only in second term). 

The discontinuity across a second order pole is derived in Appendix B. We now write down the various discontinuities 
as enumerated above: 



Disc[W^iA)l = i-2ni)8eV J 'J^^^^^^SiE - 2k) 



i2nfq 

PAfS[l/2-h{k)]' fcW5[l/2-n(fc)] 



[k - S2Eh-q - siq] [k - s^Eh-q - siq]^ 

1 2kAfS[l/2-h{k)] _ 1 k'^jAfS)' [1/2 - n{k)] 

2[k- S2Ek-q - Siq] 2 [k- S2-E^fc-g - siq] 
l fcW>S[l/2-n(fc)]' _ 1 fcW5[l/2 - n{k)] [-1 + sz^^gl 
2[k - S2Ek-q - Siq] 2 [k - S2Ek-q - siq]^ 



(5.1) 



Where, the prime denotes derivation only with respect to k. The symbol M stands for the factor [(si + S2)/2 • 
Sin{Ef^_q) + S2n{q)], while the factor (A: — q)^ ■ fc+ is represented by the symbol S. 



DisclW^iC)]. = (2.i)8eV / ^^^^|!|j^,(E - 2k) 



^ ll 2kMS[l/2-fi{k)] , 1 fc2(A/'5)'[l/2-n(fc)] 



2[E - S2Ek-q - siq - k] 2[E - S2Ek-q - siq - k] 

1 eNS[l/2 - h{k)]' 1 eMS[l/2 - n{k)] [1 + S2E'^_^] \ 

2 [E - S2Ek-q - siq - k] 2 [E - S2Ek-q - siq - kf J 

The two terms above are the result of the discontinuities at p° ^ -E = 2k. In the following we shall enumerate 
those terms that result as we take the discontinuities a.t p^ ^ E = sk + siq + S2Ek-q. 

d^kd^q [{si + S2)/2 - sin{Ek-q) + S2n(g)][l/2 - n{k)] 



/rjihr 



)6 q[s2Ek-q + Siq - sk] 
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ssik - q) ■ ks 



>5{E - Sk - Siq - S2Ek-q) 



^ s[s2Ek-q + siq - sa/c] S2Ek-q + Siq + sk 



(5.3) 



Recall that even though not explicitly mentioned there is an implied summation over all sign factors. We may now 
perform the sum over S3 = ±s to get 



Disc[\l't{C)]c = {-2ni)^e^g 



2^2 / d ^kd^q [[si + S2)/2 - Sih{Ek-q) + S2n{q)][l/2 - n{k)] 
;27r)6 q[s2Ek-q + s^q - sk^ 



X (fc - q)s2 ■ ksS{E - sk- siq - S2Ek-q) 
Finally, the discontinuity in parts D of the second self energy is given as 



Disc[Tl^|;^{D, 2 + 3)]^ = (-27ri)8eV 

sz{k - q)^ ■ ks3 

[sEk-q - sq~ S3k][sEk-q - sq~ Sik] 

S5{k - q)^ ■ ks. 



d^kd\ [1/2 - n{Ek-q - q)][n{Ek-q) + n{q)] 
(27r)6 q 

5s^,+5s,+5{E -k~ Ek-q + q) + 6s^,+6s,-S{E -k + Ek-q - q) 



[{Ek-q -qf~k 
We may now sum over S3 = ±1 to get 



6s,-5s,+5{E -k- Ek-q + q) + Ss,,-Ss,-S{E -k + Ek-q - q) 



Disc[W{D, 2 + 3)]c = (-27ri)8eV 



(k — q) ■ k+ 
X < ^ ^5{E -k + Ek-q - q) 



2 2 [ d'^kd^q [1/2 - n{Ek-q - q)][n{Ek-q) + n{q)] 



[k + Ek-q-qY 



(27r)6 q 

[k + q- Ek-q^ 



(5.4) 



(5.5) 



(5.6) 



Disc[n;;(D,l)]c = (-27ri)8eV 
f S3(fc - q)^ ■ ks 



[sEk-q + sq - szk][sEk-q + sq - s^k] 
S5{k -q)s- ks. 



d^kd^q [1/2 - h{Ek-q + q)][l - 7~i{Ek-q) + n{q)] 
(27r)6 q 



Ss^,+Ss,+S{E - fc - Ek-q -q) + 6s^,-Ss,+S{E + k- Ek-q - q) 



[{Ek-q + qf - 

We may sum over S3 = ±1 to obtain 

Disc[W{D, l)]c = (-27rj)8eV 



5s,,-5s,+5{E - fc - Ek-q -q)+ 5s,,+5s^+5{E + k- Ek-q - q) 



d'^kd^ [1/2 - hiEk-q + g)][l - niEk-q) + n{q)] 



(27r) 



(fc - q)^ ■ k+ 
[Ek-q + q- fc]' 



;diE - fc - Ek-q - q) 



(fc - q)_^ ■ k- 
[Ek-q + q + k] 



;S{E + fc - Ek-q - q) 



(5.7) 



(5.8) 



VI. PHYSICAL INTERPRETATION: TREE-LIKE CUTS 

We now begin the process of combining terms from the discontinuities of the two self energies to obtain the square 
of amplitudes of physical processes. Essentially we shall follow the method outlined by H. A. Weldon 0|. Our method 
is essentially a three step process: 
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i) collect together terms that have the same energy conserving delta functions. 

ii) reorganize the thermal distribution functions to express them as a difference of the thermal weights for particle 
emission and absorption. 

iii) reorganize the remaining momentum dependent part as the square of the amplitude of the process hinted at by 
the previous two steps. 

For easy identification wc indicate the contribution from the first self-energy topology by and from the second 
topology by 11^. We begi n with t he discontinuities where no loops are left i n the final result. These are the dis- 
continuities given by Eqs. (3J,3/7) for the first self-energy topology, and Eqs. (5^,5_^ 5_^) for the second self-energy 



topology. These discontinuities will result in physical amplitudes for three kinds of processes: photon decay, Compton 
scattering and pair creation. 



(27r)6g 

{[1 - fi{k)][l + n{q)][l - h{Eq^k)] - h{k)n{q)fi{Eg^k)} 

S3S4[k+ ■ [q - fc)+][fcs4 ■ (g - fc)s3] 
[k + q + s^Eq^kWq + Eq^k + Sik]'"'" " 



We may combine the coefficients of the same delta function from the second self-energy to get 

(Pkd^q 



Disc\r(^]{E -k-q- Ek-q) = 2 X {-2Tri)8e'^ g'^ J 



(27r)6g 

X {[1 - n(fc)][l + n{q)][l - fi{Eq^k)] - h{k)n{q)fi{Eq^k)} 
{k - q)^ ■ k+ 



[Ek-q + q-k] 



(6.1) 



A. Photon Decay and Formation. 

We begin by analyzing the terms which containing the delta function 6{E — k — q — Ek-q). The contributions to 
this from 11^ are 

-— g^[l/2 - fi{k)]S{E -k- Eq-k - q) 

S3S4[k+ ■ {q-k)^][ks^ ■ {q ~ k) + n{Eq-k) 

[k + q + ssEq-k] [q + Eq-k + s^k] 

and 

-— /[1/2 - fi{Eq_k + qWE -k- Eq_u - q) 

S3S4[k+ ■ {q-k)^][ks^ ■ (q - k)^J[l - n{Eq_k) + n{q)] 
^ [k + q + ssEq-k] [q + Eq-k + Sik] ' 

Note that 

n{Eq-k +q) I - n{Eq-k) + n{q) = n{Eq-k)n{q) 
using the above identity we may combine the two terms and rewrite the distribution functions to give, 

Disc[Ii\E -k- Eq_k -q) = {-27n)8e^g^ J '^'^'^''^ 



(6.2) 



S{E-k-Eq-k-q)- (6.3) 



;S{E-k-Ek-q-q). (6.4) 



The overall factor of 2 is the ratio of the symmetry factor of this diagram to the denominator obtained from 
perturbation theory. Note that we obtain the same form of the distribution functions, this indicates the generic 
structure of heavy photon decay and reformation. In the distribution function factor, terms like 1 -I- n{q) indicate 
Bose-Einstein enhancement in the emission of a gluon. The 1 is from spontaneous emission, and the n{q) represents 
stimulated emission of a boson into a thermal bath. Terms like 1 — h{k) represents the "Pauli blocked" emission of 
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a quark of momentum k into the thermal bath. The product of the three factors [1 — + n((7)][l — n{Eq-k)] 

(along with the phase space integral and delta function) can thus be interpreted as the statistical factor associated 
with a heavy photon outside a thermal bath decaying by emitting a quark, antiquark, and a gluon into a thermal 
bath. Subtracted from this is the factor n{k)n{q)n(Eq-k)', this represents the formation of a heavy photon from a 
quark, antiquark and a gluon, all three emitted from the thermal bath (Fig. ^). The photon subsequently escapes 
from the bath without further interaction. 




FIG. 3. Heavy photon decay and formation. 



To convert the above expressions into cross sections for heavy photon decay and reformation, we start by first 
defining a new four- vector w = {w,w), such that w — E ~ k ~ q{m order to avoid confusion we introduce the notation 
of four- vectors as bold face characters). This relation is indicated by the one dimensional delta function. To obtain 
the probability of photon decay, we need to generalize the delta function to a four-delta function. We thus need to 
generalize the definition of w: 

w = p k q. 

Where, p — (£',0,0,0) is the mass of the off-shell photon. As denoted by Fig. ^. k, q and w are all on shell. The 
above relation also implies w — — fc — q. We may set h'^ ~ k = \k\ and q!^ = q = \q\. Now, requiring that w be on 
shell imposes the condition that 

[E - k - qf ^ k^ + q^ + 2kq cos 9 



E{E -2k- 2q) = -k • q. 



Where 9 is the angle between the three vectors q and k. Using the above relations we many now rewrite the 
discontinuity obtained from iP. In the numerator of the integrand we notice the factor {k — q)^ ■ k^, this may be 
changed appropriately by setting k ^ —k in the integrand. Noting that (— fc)^ = —k-s we get the above factor as 
—111+ ■ k-. Introducing the standard denominators 2k, 2w and factors of 27r we obtain Eq.(6^) as 

D^.sc[U\E -k-q- £,_,) = -2^ x 8e'g' J ^^^^^2^)H\p - k - q - w) 



{2T:f2q2k2w 
I [1 — n(fc)] [1 + n{q)\ [1 — n{w)] — n{k)n{q)n{w) | 



q-kf 

We now split the above integrand into two parts and in one of them switch w ^ k. Note that 
wk + w ■ k — {E — 2k){E — 2w)/2, finally gives the above discontinuity as 

(Pkd^qdPw 



(6.5) 

-w_|- • k_ = 



Disc[W]{E -k-q 



{2TTf2q2k2w 



(27r)'^(5*(p-k-q-w) 
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X |[1 — r;,(fc)][l + n(<z)][l — n{w)] — h{k)n{q)n{w)^ 
X 32e2g2 



E-2k E-2w 
E-2w ^ E-2k 



(6.6) 



We now perform the same procedure on the corresponding discontinuity from H^, to get 



Disc[n^]{E -k-q- Ek-q) = -i x 8eV 



(27r)92q2fc2u) 



8(2^)4r(p-k-q-- 



X |[1 — n(fc)][l + ?^((7)][1 — h{w)] ~ n(fc)n(q)n(w)| 
[/c + g + S'iw\ [q + w + S4fc] 



(6.7) 



The part of the integrand besides the distribution function part (depends on the angle between k and g, will be denoted 
as the matrix part) may be expanded by summing over S3, S4 as 



k • w 

kw 



k+ • w+ 



k+ • w 



k_ • w+ 



k_ • w_ 



{k + q + w)"^ {k + q)^ — (q + w)^ — k'^ [q^ — (k — wy 



Using the relations 



= k 



f+ = -~{1/2)[E -2k][E -2w], 



k+ • w+ = k_ • w_ = {E/2)[E - 2q], 

and the relation imposed by the delta function {i.e.,E = k + q + w) we can simplify the matrix part to give 

E{E - 2q) 
[E-2w][E-2ky 

substituting the above into the expression for Il^and then combining the results from Il^and Il'^we get 

d^kd^qd'^w 



Disc[Il\{E -k-q-E, 



k-qj 



{2TTf2q2k2w 



(27r)'*(5*(p-k-q-w) 



I [1 — n{k)] [1 + n{q)] [1 — n{w)] — h{k)n{q)h{w) | 

2 2 £'-2fc E-2w 2 E{E~2q) 
E-2w^ E-2k ^ [E - 2w] [E - 2k] 



X 32e^g 



(6.8) 




FIG. 4. Heavy photon decay at first order in a and as- 



Photon decay into a quark, antiquark and a gluon at first order in the E.M. and strong coupling constant can occur 
by two types of Feynman diagrams p^ ] as shown in Fig. ^ The matrix element for the first diagram may be written 
as Ml — Miei_i{p), where 
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for the second diagram 



= u{k)igYe*p{q)''-j^-^ieYv{w) 



Taking the product A4^^A^i^ and summing over the spins and colors of the quark, antiquark and the gluon gives 



M^Mi^ = -32eV 



E-2w 
E-2k' 



Similarly 



bj — 2w 



(6.9) 



(6.10) 



Notice that as the three 3-vectors fc, g, w form a triangle, E — 2k — w + q — k\s always positive. By the same 
argument E — 2w and E — 2q are also positive. We thus note that A^|''A^i^, A^2''-^2;j are negative. This is to be 
expected as the square of the full matrix element is positive, where from the sum over the photon's spin we get 



The cross term is 



2 2 2E{E - 2q) 
-6Ze g 



[E - 2w] [E - 2k] ' 



(6.11) 



Comparing the above three equations with the result from the loop calculation (Eq.6^, 6^) gives us the relations: 

ffkd^qd^w 



Disc[n";]{E-k-q-Ek-q) = 



(27r)92g2fc2w 



(27r)M^(p-k-q-w) 



|[1 - n{k)][l + n{q)\[l - n{w)] - n{k)n{q)n{w)^ 



(6.12) 



Disc[n'';:]{E - k - q - Ek-g) 



d^kd^qd^w 
{2nf2q2k2w 



(27r)454(p-k-q- w) 



|[1 - n(fc)][l + n{q)][l - fi{w)] - n(fc)n(g)n(w)| 



and hence we get the relation written down by Weldon 

d'^kd'^qd'^w 



(27r)92(j2fc2w 

X |[1 — n(fc)][l + n((7)][l — n{w)] — n{k)n{q)n{w)^ 



(6.13) 



(6.14) 



where M. = M^ef_i{p) — Mi + M2 is the full matrix element of heavy photon decay. 
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B. Compton Scattering. 



The an alysis for Co mpton scattering is slightly more tricky. Note that there are two sets of terms from Eqs. (3.6,^) 



and Eqs. ([3.4,5.6,5 



that may lead to Compton scattering. One appears with the delta function d{E + k — q — Ek-q) 
and the other with the delta function 6{E + Ek-q — k — q). The delta functions can be turned into one another simply 
by replacing k ^ k + q, followed by g ^ —q. One notes on performing this operation that the rest of the integrand 
looks rather different. This happens as there are 4 topologically distinct diagrams that may fall under the category 
of Compton scattering ( it is well known that that for a given in state there are two diagrams that lead to Compton 
scattering; there are four here as we sum over the possibilities of the incoming fermion being a quark or antiquark). 
Let us consider the contribution from H^: 



Disc[U^]{E + k- Eq_k - (?) = {2Tri)8e^g^ 

X |n(fc)[l 



(Pk(Pq 
"(2^ 



fi{Eq-k)][l + n{q)] - [1 - fi{k)]fi{Eq-k)n{q)^ 
{f^k)^][k,-{q^k)^^ 



[-k + q + sj,Eq^k] [q + Eq^k + s^k] 



5{E + k 



<!)■ 



(6.15) 



For the contribution from H^, recall that we have an overall factor of two on each of the results of Eqs. (5.4,5^, 5.8 ) 
coming from the overall symmetry factor of H^being double that of H^. We take half of the contribution from the 
-q) term, and half from the 5{E-\'Ek^q — k — q) term, and in the second contribution change /c k + q, 



S{E + k~Ek-q- 

followed hy q — 



-q. This gives the total contribution from 11 as 



Disc[U'^]{E + k- Eq^k -q) = {27Ti)8e^g" 



d^kd^q 
{2TT)^q 



{n(fc)[l - fi{Eq^k)][l + n{q)] - [1 - n(fc)]n(i;,_fc)n(g) } 
k--{q- fc)+ ■ {q- k)^ 



[k- 



E, 



'q-k\ 



[k-q + Eq.k? 



SiE + k-Eq.k-q). 



(6.16) 



Notice that the combination of distribution functions appearing in the curly brackets are identical. The product 
of the three factors n{k)[l — ri(i?q„fc)][l + n{q)] has the interpretation of an incoming quark(or an antiquark) from 
the medium fusing with the photon coming in from outside the bath, resulting in an gluon and a quark(or antiquark) 
going into the medium. Subtracted from this is the product [1 — h{k)]h(Eq^k)'>T-{g) , which has the interpretation 
of an incoming quark(antiquark) from the medium fusing with an incoming gluon from the medium, resulting in a 
quark(antiquark) going back into the medium, and a virtual photon which leaves the medium (Fig. |^). 




FIG. 5. Quark Compton scattering. 



To convert the above expressions into cross sections for Compton scattering, we again define the new four-vector 
w = {w, w), such that w = E + k — q. This is generalized to 

w = p + k q. 

as a result w = k — q. Now, requiring that w be on shell imposes the condition that 
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{E + k - qf = + - 2kq cos 9 



E{E + 2fc - 2(7) k • q. 



Using the above relations we many now rewrite the discontinuity obtained from 11^. In the numerators of the 
integrand we notice the factor (k — • fc_, which may be written as w+ • k_ = (£' + 2k){E — 2w)/2. We introduce 
the standard denominators 2k, 2w and factors of 2tt and perform a similar set of operations as for photon decay to 
obtain the full result for Compton scattering as 



Disc[U]{E + k - q ~ E, 



k-q) 



(PkcPqd^w 



(27r)92q2fc2w 
|n(fc)[l + n(g)][l - n(ti;)] - 

E + 2k E-2w 



(27r)4^'*(p + k-q- w) 



2„2 



X 32e^5 



i(fc)]n(q)n(?i;)| 
E{E - 2q) 



E-2w E + 2k [E~2w][E + 2k] 



(6.17) 







FIG. 6. Compton scattering at first order in a and as. 



Recall that in H^there is another term, the coefficient of the delta function 6{E + Eg-k — k — q), that leads to 
Compton scattering. Also, in the Compton scattering contributions from iP, we only used a half of both the terms. 
Following almost the same method as above, one can demonstrate that the form of the contribution from these terms 
is almost the same as above with k and w interchanged. In it, one may interchange ix; — > to get the same contribution 
as Eq. (6.17); hence doubling the total contribution from Compton scattering. 

Compton scattering by an incoming photon of a thermal medium of quarks and antiquarks, at first order in the 
E.M. and strong coupling constant can occur as a result of four processes as shown in Fig. ^. The matrix element for 
the diagrams may be written as Ain = Ail^e^{p), where 



(p + k)^ 



■167^6,, (p)ii(k) 



(w - p)2 



'e;(q)u(k). 
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The amplitude for the third and fourth diagram can be obtained from the two amphtudes above simply changing 
u ^ V. Taking the products and summing over spins and colors (remember diagrams 1 and 2 interfere with each 
other, and so do 3 and 4), we get 



M*/Mi^ = 32eV 



, E-2w 
E + 2k' 



(6.18) 



E-2w 



(6.19) 



Once again, note that A4l^\Mi^,Ai2'^M2n are negative. This is because E — 2w — q — k — w is always negative 



due to the triangle condition mentioned in the previous subsection. The cross term is 



<*A'w _oo„2„2 2E{E-2q) 



[E - 2w] [E + 2k] ' 



(6.20) 



Comparing the above three equations with the result from the loop calculation (Eq. (6.17)) gives us the relation: 

mscW,]iE + k^q- E,.,) . . / ^^^i2^mp + k - q - w) 

X |n(fc)[l + n{q)][l - h{w)] - [1 - h{k)]n{q)n{w)^ 



(6.21) 



Once again we note the interesting fact that in this gauge the mixed terms A^2^A^i^ + A^i^A^2/j are always given 
by n^and the square terms Ml^Mi^^ + M2^M2f_i are furnished by H^. 



C. Pair Creation 







FIG. 7. Pair creation at first order in a and as. 



The analysis for pair creation(often referred to as photon-gluon fusion) is almost identical to the two previous 
sections. I ts contribu tion is furnished by the only remaining delta functions in Eqs. (3.6|,3.7) in the first self-energy, 
and Eqs. (5^,^,5^) in the second self-energy, i.e., S{E + q — k — Ek-q). We simply state the results here: pair 
creation can occur through two types of processes, and has the discontinuity in the total self energy as 



Disc[\l\E + q - k ~ Ek- 



d^kd^qd^w 



(27r)'*r(p-k + q-w) 



(27r)92q2fc2w' 
X |[1 - n{k)]n{q)[l - n{w)] - [1 + n{q)]n{k)n{w)^ 



X 326^5^ 



E-2k E-2w 



E{E + 2q) 



E-2w E-2k [E- 2w] [E - 2k] 



(6.22) 
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VII. PHYSICAL INTERPRETATION: LOOP-CONTAINING CUTS 



We now analyze the various discontinu ities of Il^and H^which contain loops. We start with the discontinuities of 
n^. These are given by Eqs. (3.3 3.4,pT5). We note that there are two terms with the delta function S{E — 2Ek-q) 
these correspond to the cut of Fig. |^. There is, also, one term with the cut 5{E — 2k), this corresponds to the cut of 
Fig. |. 




FIG. 8. Photon decay at one loop corresponding to the cut 5{E ~ 2Ek-g). 



One may be satisfied with this interpretation of the cut diagrams and not proceed further. A recent paper [[l2| , 
however, has drawn attention to the fact that one can obtain a somewhat different interpretation of these diagrams, 
in terms of interference between simple tree like diagrams and diagrams containing particles called "spectators". 
Spectators are essentially on-shell particles from the heat bath that enter with the in-state and leave with the out- 
state without having interacted with the the rest of the "participants" . 




FIG. 9. Photon decay at one loop corresponding to the cut 5{E ~ k). 



We start by summing over the variable si in Eq. 3.3. This immediately gives two terms, distinguished by the 
combination of distribution functions they carry: 



2„2 



Disc[2U'^{A)]ai = (-27rij»e-5 

X S2 



d^kd^q[k+-{q^k)J[k^-{q-k)^ 
(2^)6 q{E - (S2 - s:i)Eq_k) 
l-2n(fc)][l-2n(£;,_fe)]g 



(fc -1- s:iEk- 



-5{E-2k). 



and 



Disc[2W{A)]a2 = (-27ri)8eV 



d^kd\[k+-{q'k)J[k^.{q-k)J 



(27r)6 

X (-S2)[l -2n(fc)][l + 2n(g)] 



q{E - (s2 - S'i)Eq^k) 

Sjk + Eq^k 

{k + s^Ek-qY ~ q^' 



6{E - 2k) 



In Eq.7.1, if we replace S2 



-S4,S3 — s- — s, followed by fc q^ k we get 

d^kd'q [{q'^k)_ ■ kMo^k)^ ■ ks] 



Disc[2Il^^iA)U = {+27n)8e'g 



{2n) 



q{E - {s - S4)Eq^k) 



[l-2h{k)][l-2fi{Eq.u)]q ,,^ „P . 
" {Eq.,-sk)--q- - 



(7.1) 



(7.2) 



(7.3) 
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The above is exactly the same as Eq.(3.4). This is to be expected as the two cuts should in principle represent the 
same diagram up to a shift in momenta. We thus double this contribution and focus on it. It represents photon decay 
into two quarks with quark emission and absorption from the final state quarks. The other part from Eq. (7.2) along 
with Eq. will represent photon decay with giuon emission and absorption off the external quarks. 



A. Photon decay with quark emission-absorption off vertex and final state. 

We begin by summing over the remaining sign variables S2, S3 in Disc[2Il'^{A)]ai to get 

/d^k(fia 
——^[1 - 2n{k)][l - 2n{Eq^k)]5{E - 2k) 
(2tt)'=' 

[k+ ■ iq^k)^][k- ■ jq^k)^] [k+ ■ {q^k)_][k. ■ jq^k)^] 
E[{k + Ek-,Y - g2] [E + 2E,^k)[{k + E^-^Y ~ q^] 

[k+ ■ {q^k)A[k- ■ {q^k)_] [k+ ■ {q^k)_][k- ■ {q^k)_ 



{E - 2Eq^k)[{k - Ek-,? - q^] E[(k - Ek-qf - q^] 



(7.4) 



As in the previous section, the distribution functions will be reorganized to allow for an interpretation in terms of 
thermal weights for particle emission and absorption. In the first two terms we define the new lightlike four-vector w 
such that w — q — k. In the last two terms we define w such that w = —q + k. This allows us to change the variable 
of integration as (Pq — > (Pw, as fc is a constant as far as the q integration is concerned. We may also redefine the 
distribution functions as 

[1 - 2n(A:)][l - 2n(w)] = [[1 - n(fc)][l - fi{k)] - n(A:)n(fc)j [[1 - n(w)] - [n{w)] 

The first set of factors in the larger square brackets has the usual interpretation |^ of the thermal factors that are 
associated with the probability of particle emission into a heat bath or particle absorbtion from a heat bath. In this 
case they carry the obvious meaning of: [emit fermion of energy k] [emit fermion of energy k] — [absorb fermion of 
energy k] [absorb fermion of energy k] . The reader will note that unlike the self-energy cuts considered in ^ or those 
of the previous section, the two cut diagrams that will result from this imaginary part of the self-energy will not be 
symmetric, in the sence that it will be the interference between a diagram with a loop and a simple tree diagram. 
The thermal factors discussed above will be the same for either diagram as they pertain to the quark and anti-quark 
that emanate from the decay of the photon (or those that combine to form the photon) . Both amplitudes that result 
from this imaginary part contain this process and thus have identical thermal factors. 

The second set of thermal factors has a new interpretation. These thermal factors pertain to the particles in the 
remaining loop and thus are germane to only one of the two interfering amplitudes. We will demonstrate that these 
signal the difference of two amplitudes: that for the emission of a quark or anti-quark into the bath and its subsequent 
absorption from the bath, and vice-versa. Thus the second set of distribution functions is to be understood as: [emit 
fermion of four-momentum w] [absorb the same fermion of four-momentum w] — [absorb fermion of four-momentum 
w] [emit the same fermion of four- momentum w] . 

The process of emission of a fermion of four-momentum w into a bath followed by its reabsorption is formally 
achieved by the action of creation and annihilation operators on the bath state [riw) i-e, 



J^\n^). (7.5) 



The reverse process, i.e., the absorption of a fermion from the bath and it subsequent re-emission into the bath is 
formally achieved by the action of annihilation and creation operators on the bath state i.e., 

a^a\h^) = (n^^\n^). (7.6) 

The discontinuity of the self-energy will represent the amplitude of a particular process multiplied with the complex 
conjugate of another. In one of these processes the above mentioned fermion will perform the emission and absorption 
procedure referred to above. In the other amplitude, as we will show shortly, it will simply enter and leave without 
having interacted with the rest of the particles. Due to this reason, it has been referred to previously (see |l^,^) as 
a spectator. 
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FIG. 10. Interference between diagrams of different order in as. The diagrams on the left indicate 2 — » 3 reactions lilce 
7(5 — > qqQ (where the Q indicates that the incoming and outgoing quarks are identical). The diagrams on the right indicate 
the complex conjugate of Born term photon decay with a comoving quark spectator i.e., (7 — > qq) ® {Q ^ Q). 



We introduce the usual denominators of 2w, 2k to get 



Disc[m]4 = 2 X (-27ri)8e2.g2 



-8{E - 2k) 



+ 



{2Tr)^k'wkw 
[1 - n(k)][l - n{k)] - n{k)h{k) 

[1 — n{w)] — [n{w)] 

[k_|_ • w+] [k_ • w+] [k+ • w_] [k_ • w+] 

E[{k + w)"^ - q^] 

[k+ •w_][k_ • w+] 



{E + 2w)[{k + w)"^ - q"^] 
[k+ • w+][k_ • w+] 



{E - 2w) [{k - w)2 - g2] E[{k - m;)2 - g2] 



(7.7) 



We now introduce the new four-vector k^ = (fc, —k), and generalize the delta function to a four delta function. We 
then combine the first two terms and the last two terms to write 



^ ''^ ^ J {2wf8kwkt 



16(27r)''(5'*(p + w - k - kb - w) 



+ 



[1 - h{k)][l - n{kb)] - h{k)h{kb) 

[1 — h{w)] — [n{w)] 

[kfe ■w][ka ■ (p + w)] 
[{E + w)'^ - w^][{k + - 

[kg ■ w] [kb ■ (w - p)] 

[(w-i;)2 -w2][(fc-«;)2 -g2] 



(7.i 
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The above, has the interpretation of Fig. |0[ This indicates the interference between two diagrams of different order 
in couphng constants. Let the matrix elements of the two tree-level diagrams with two propagators be denoted as 
A4i = A^^e^(p) and A^2 = -^2^m(p)- matrix element of the term in brackets is simply denoted as m^e^(p). 
Where the dotted line called the spectator is simply a product of Dirac delta functions over four momenta and 
Kronecker delta functions over the spins and colours of the incoming and outgoing fcrmions denoted by Wa and Wh 
(here, for brevity we indicate all the different quantum numbers, both continuous and discreet, of the incoming and 
outgoing particles by a single label). 

It is now simple to verify that the result obtained in Eq. (^^) can be written as 



Discm 



d^kd^wd^kb 



- n{k)n{ki,) 



{2TTf%kwkb 

[I ~ h{k)][l - n{kk)] 



k k. 



{1 — n(ti;)} — h{w) 



(7.9) 



Where the Kronecker and Dirac delta functions over the fermions Wa and Wh have been used to set Wa = Wb = w. 
The factor of 2 preceeding the interference matrix elements is due to the fact that a similar process may be obtained 
by replacing an incomming quark spectator with an anti-quark spectator. 



B. Photon decay with gluon emission-absorption from final state quarks. 



This term receives contributions from Disc[2'n'^{A)]a2 and Disc[2Il'^{D)]b. The fate of this discontinuity is essen- 
tially similar to the previous section and results once again in the interference of tree level diagrams of different order. 
There are two sets of diagrams with two propagators here as well, the difference being that the incoming, outgoing 
particle with the same set of quantum numbers or in other words the spectator is a gluon. We once again introduce 
the on-shell four-vector w, such that w = k — q. We use this to change the variable of integration in Disc[2Il'^{D)]i,. 

Disc[2Il^(A)]a2 we relabel the dummy variable k ^ w. Both discontinuities give essentially the same contribution 



thus the total discontinuity from such processes (Fig. 11) is given as 





"b 

FIG. 11. Photon decay with spectator gluon (the Born term with spectator gluon is implied, see Fig. [lo| ). 



d^wd^qd^Wb , 



{2iTfSkwkh 
1 — fi(w)][l — n[wb)] — h(w)fi(wb) 

X [l + n{q)] + [n{q)] 



32(27r)''r(p + q - v^r - v^Tfe - q) 
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[wfc ■ (q - w)][Wa ■ (q - Wfc)] 

[(g-w)2 -fc2][(g + w)2 -fc2 



(7.10) 



Which is once again equal to 



^^^^P^^l^ = ' J Wnk^^ <5 (P + q - w - w, - q) 

[1 — — fi{wb)] — n{w)n{wb) 

[1 + n{q)] + [niq)] 



(7.11) 



Where m represents the same process as in the previous subsection. The amphtudes A4i and A42 represent the 
processes of Fig. |ll|. The interpretation of the first set of distribution functions is the same as before i.e., emission and 
absorption of two particles of energy k. The second term has the interpretation of a gluon spectator exactly identical 
to that of the quark spectator in the earlier subsection, but with Pauli factors replaced with Bose factors. Note that, 
unlike in the previous section, there is no factor of 2 preceeding the matrix elements as the spectators are gluons. 



C. Photon decay with quark and gluon emission-absorption off the same quark line 



We , now, begin the analysis of the last loop-containing cut. This is essentially given by the discontinuities of 
Eqs. (13,13). Comb ining these two discontinuities, and writing k = E — k in the denominators of the terms from 
£'zsc[n^(v4)]a (note that we have to double this contribution as it emanates from the second self-energy diagram 
which has a symmetry factor of 2 more than the first self-energy diagram), we get 

Z^z^cP^e = i2.^)SeV J '-M^^^^^SiE - 2k)[l - 2n(fc)] 
f 2kJ\fS k^{AfSy 



][E~ S2Ek-q ~ Siq - fc] [E - S2Ek-q - Siq - k] 

k^^S[l + S2E',_^] \ 
[E^S2Ek-q-siq-k]^j- 

(2..)8eV / ''7/;;;'''^ ^(i^-2fc)[l-2n(fc)] 



d_i k^AfS 



x^T n^^-TT— ^ ■ (7.12) 



dk 



\^[E - S2Ek-q - Siq - k] 



The above term does not readily admit a physical interpretation, however, the infrared limit will be evaluated with 
the above expression as the starting point as it is formally correct. To try and obtain a physical interpretation from 
the expression given above, an integration by parts is performed to obtain the discontinuity as 

mscK], = (2..)8eV / - - 2HE/2)] 



[E - S2Ek-q - Siq - k] 



{2m)%e^g^ 



d^kd^q MS 



(27r)6g [k - S2Ek-q - Siq] 
X ^^2S'{E - 2k)[l - 2n{E/2)]y (7.13) 

Where, we have used the property that 
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dSiE - 2k) ^ ^d5[E^ 2k) ^ _^dSiE - 2k) ^ ^^^^ 



dk 



d{2k) 



d{E) 



Interestingly, as an aside, we note that one may still obtain a physical interpretation of the above term in terms of 
spectators with retarded propagators. To obtain this, we expand the factor -siq ] summing over si and 

S2- Here, as expected, we will obtain a part dependent on Bose distribution functions and a part dependent on Fermi 
distribution functions. We will illustrate the physical interpretation using the part containing the Bose distribution 
functions. We begin by writing the delta function in Eq. (7.13) using the following regulator: 



5{x) = lim 



In this regulation scheme, we obtain 



5'{x) 



5{x) 



X + ie 



+ iTr6'^{x) 



Substituting the above relation in Eq. (7.13) we get 

d^kd^ 



^^(-2)<5(£;-2fc)[l~2n(i?/2)] 

+ iTT5{E-2k) 



NS 



[E - S2Ek-q - siq -k]y{E- 2k) 
We now write unity in the form of an integral as 

1^ I dco"i{5(Lj" + fc) + 5(c^o + fc)} 

duj^kdioj^"^ ~ k^). 

Substituting Eq. ( |7.17| ) in Eq. ( [7.16D , we obtain 

Disc[W;]e = (27ri)16eV J dio''{E/2)S{Lo'^^ - {E/2)^) [ ^^^i-2)S{E - 2k)[l - 2h{E/2)] 
JVS f 1 



[E - S2Ek-q - siq -k]\{E- 2k) + le 
- 27ri(n(tj°) - 1/2)5{E - 2k) \ 



(27r)6g 

+ 2TTin{uj^)5{E - 2k) 



d^kd^q 

(277)6^ 



(27ri)(-16eV) j dio^ 

EJVS I 1 

[E - S2Ek-q - Siq -k]]^{E- 2k) + ie 

- 2ni{n{uj^>) - e{~J>))S{E -2k)\. 



(5(cj° - {E/2)^)SiE - 2k)[l - 2h{E/2)] 



+ 2TTin{uj°)S{E - 2k) 



(7.14) 



(7.15) 



(7.16) 



(7.17) 



(7.18) 



In the above we have simply added and subtracted the factor 2TTih{uj'^)S{E — 2k) inside the curly brackets and 
rewritten 1/2 as 6'(— w"), as the rest of the integrand is an even function of 
We now introduce the three vector part of uj as 



(5(w°^ - {E/2)^)SiE ~k^k)= I d^uj5{uo^'^ - {E/2)^)5{E -k~ k)S^{~k ~ uj) 



d^ujSiJ^'^ - \Q\^)8{E - fc - M|)<53(-fc - w). 



(7.19) 
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As stated previously we now concentrate on the part of MS /\k—S2Ek-q—siq\ which depends on the Bose distribution 
function. This gives 



d^UJ 



2^5(w*' - \uj\^){2TTf5^{i> - k - w) 



(27r)4 J (27r)62g2fc' 
{1 - n(A:")}{l - - [{1 + n{q)} + n{q) 

k{E~k~q)-k-{k-q) k{E-k + q)-k-{k-q) 



X (-64e^g2)£;2 



d^kd^q 



27r5(t^°^ - Pn(27r)'*(5''(p - k - w) 



(27r)4 7 (27r)62g2fc 
{1 - n(fc")}{l - - n(fc°)n(w")l f{l + + n(q) 



(7.20) 





FIG. 12. Photon decay into a pair. The quark then emits or absorbs a quark or a gluon. 



Where i i/)Aj^{uj) is the retarded propagator. One may note that the integrand in the above equation is simply the 
interference matrix elements of the first {A4i — e^(p)Ali^p) and second diagram (A^2 = e''(p)A^2,/j)of Fig. |lj and 
the Born term (m — e^(p)m^') with a gluon spectator. A similar interpretation may be obtained for the the third 
and fourth diagrams of Fig. n3 in terms quark spectators. Howev er, a. s the above equation is not mathematically well 
defined; it will not be used in evaluating the infrared limit. Eq. ( 7.1^ ) will be used instead. 
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VIII. INFRARED BEHAVIOUR 



We now examine closely the infra-red and coUinear singularity structure of the terms enumerated in the two sections 
above. We will examine the infra-red behaviour in the limit of heavy dilepton production from a plasma of massless 
quarks i.e., E » T. There are essentially five terms: 

F: Photon Gluon production denotes the reaction q + q ^ g + J 

C: Compton like reaction between a gluon and quark/anti-quark g + q q + J 

D: denotes the three body fusion to form the photon g + q + q ^ J 

V: denotes photon formation from vertex corrected quark, antiquark. 

S: denotes photon formation from self-energy corrected quark, anti-quark. 

The full imaginary part of the two loop self-energy may be schematically written as 

2Imn|l = J dw[n{w)[Fiw) + Dg{w) + Vg{w) + Sg{w)] 

+ n{w)[C{w) + Dg{w) + Vq{w) + Sq{w)]y (8.1) 

The first four terms represent the part of the terms mentioned above which are proportional to the gluon distribution 
function. The last four terms are those proportional to the quark/antiquark distribution function. This is essentially 
the same notation as used by the authors of H. We now compute these contributions in turn. 



A. Self-Energy correction Sg &: Sq 



The self-energy correction is essentially given by Eq. ( 7.12 ), i.e., 

dkd0d(p sin 0d^q 



Disc[n';]e = (2^i)8e2g2 



E 



(27r)69 



S{E ~ 2k)[l - 2n{k)] 



d 



dk [E - S2Ek-q - siq - k] ' 



We concentrate, first, on the sum S = J2si 



This may be expanded as 



5 = i 

w 



S2 [E-S2Ek-q~siq~k] ' 

[1/2 - n{w) + 1/2 + n(g)]w+ ■ k+ [1/2 - n(w) - (1/2 + n(g))]w+ • k_ 



E — k — w — q E — k + w — q 

(1/2 - n(w)) + 1/2 + n(q)]w+ ■ k+ [-(1/2 - fi{w)) - (1/2 + n{q))]xv+ ■ k_ 



E—k—w+q E—k+w+q 
We now concentrate on the terms proportional to 1/2 + n{q), i.e.. 



Sg^k 



\ + n{q) 



2k(E -k-q) -2w-k 2k{E - k + q) - 2w ■ k 



{E-k-qY- w2 



{E -k + qY - w"^ 



(8.2) 



(8.3) 



(8.4) 



Introducing the variables a = E — 2k — 2q, (3 = E — 2k + 2q and y = cos 9{ where 9 is the angle between k and q), 
we obtain 



Sg^k 



1 



n{q) 



{2k-E)a 



{2k - E)(i 



Ea + 2kq(l + y) Ef3-2kq{l-y) 



(8.5) 



Dropping the ^ ahead of the gluon distribution function we obtain the matter part of Sg . Using only this part we 



obtain (performing the unimportant angle integrations) 



8eV 
(27r)3 



dqn{q)Sg{q) 



8eV 
^2^ 



dk[l - 2h{k)]Sik - E/2) I dqq I dy—Sg,rnat 



d 



(8.6) 
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The limits of the y integration are the locations for the onset of coUinear singularities, these are shielded by removing 
a small part of phase space e i.e., the y integration is performed within the limits — 1 + e— > 1 — e. The results will 
now depend on e. This gives the result as 



8eV 
(27r)3 



dqn{q)Sg{q) 



(27r)3 



dqn{q) 



Aq — 4q log 



■7) 



In the above, the term 2n{E/2) has been dropped, as we are interested in the heavy dilepton limit where E >> T 
and as a result n{E/2) 0. Thus, we get 



Sg{w) = —Aw — Aw log 



We now concentrate on the terms, in Eq. (B.3), which are proportional to the factor 1/2 — h(w). Following a similar 
procedure as above we obtain 



2„2 



Thus, giving us the relation 



2„2 



dwh{w)Sq{w) 



dwn{% 



Aw — Aw log ( — 



Sq{w) = —Aw — 4ii;log I - 



(8.9) 



(8.10) 



Note that in both the expressions for Sq and Sg there is a log (|) term which blows up as e 0. This is a collinear 
singularity. We shall allow e to vanish only when all the different contributions to the dilepton rate have been added 
together. 



We concentr ate f irst on th e term proportional to the fermionic frequency i.e., Vq. This vertex correction is essentially 
given by Eqs. (7T, 7^, 1_A). The first two need to be doubled, as mentioned before in section 6. Extracting only the 
part proportional to the fermionic distribution function h{w), we obtain the Vq integral as 



8eV 
(27r)3 



dwn{w)Vq{w) = 2 X (-27ri)8e^.g- 



(k4 



d^kd^w 



[E — (s2 - S3)w] [{k + s^wy 



{[l-2n{k)][-Mw)]} 
S{E - 2k). 



(8.11) 



Performing the sum on S2, S3 and setting y = cos0(where 6 is the angle between k and w), we can perform one of 
the integrations with the help of one of delta function to get 



{2n 



V , 32eV 
— J / dwn[w)Vq{w) — 



(2^)3 



dww{~n{w)) 



nl-c 



-l + € 



dy\ 



E 



E 



1 + y 

E + 2w^ E -2w\\-y 



E + 2w E -2w 



■(1 + y) 



.12) 



Note, once again, that the limits of the final angular integration y signal the onset of collinear singularities. These 
are, once again, regulated by removing the small part of phase space e. At this point we introduce the condition 
that the limit of interest is for dilepton mass much larger than the temperature i.e., E >> T. The presence of the 
distribution function h{w) on the energy w severely restricts the contribution from regions where w >> T to the 
integral. Thus, the dominant contribution to the integral is from the regions where w << T or w ^ T. Hence, in the 
integral we may make the approximation that w << E and expand the factors in the square brackets to linear power 
in w/E. This finally gives 
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86^52 



(27r)3 



dwn{w)Vq{w) 



8eV 



dwh{w) 



8w + Sw log ( - 



Thus we obtain that 



Vq{w) = -8w + Swlog - 



.13) 



(8.14) 



Following almost a similar method as above we may obtain Vg from Eq. (^^) (with an overall factor of 2 as there 
is another cut which gives an identical contribution) as, 



2E 

Vn (w) — Aw log ( - 

w ' ' 



(8.15) 



Once again, note that both expressions demonstrate a coUinear divergence as e — > 0. The term Vg also displays an 
infrared divergence as w — > 0. 



C. Photon formation from quark, antiquark and gluon Dg & Dq 

The reverse reaction to this process represents heavy Photon "decay" in to a 997. Due to this reason, the process 



is denoted by the letter D [g| . The full decay contribution is given by Eq. (|6 



as, 



Disc[W]{E - k - q - E, 



k-q) 



,^ ^9. (2^)^5^(p-k-q-w) 
{2'K)^2q2k2w 

|[1 - h(k)][l + n{q)][l - h{w)] - h{k)n{q)fi{w)^ 

, E{E - 2q) 
' [E - 2w\[E - 2k] 



X 32e2g2 



E-2k E-2w 



E-2w E-2k 



(8.16) 



In the above equation, note that if three of the delta functions are used to set w — —k — g, then the remaining delta 
function imposes the condition that 



E = k 



\/k^ + g2 



2kqcos9 



As mentioned before, we work in the limit E >> T, m this case the delta function can be satisfied by the following 
regions of phase space: 

a) k ^ E, q E and hence w ^ E; in this case all the distribution fmictions n{q), fi{k), h{w) 0, and, thus, so do 
products of distribution functions. 

h) k ^ T « E, q ^ E and hence w ^ E\ in this case n[k) ^ 1. However n{q), h{w) 0, and so do products of 
distribution functions. 

c) w ^ T << E, q ^ E and hence k ^ E\ in this case n{w) ^ 1. However n{q),n{k) — » 0, and so do products of 
distribution functions. 

d) q ^ T << E, k ^ E and hence w ^ E\ in this case n{q) ^ 1. However n{w),n{k) 0, and so do products of 
distribution functions. 

Contributions from b) and c) will give us Dq, d) will give us Dg, while the contribution from a) is negligible 
in comparison. We begin by calculating Dq from the regions b) and c) of phase space. Here we can ignore all 
combinations of distribution functions containing n{q). As before, we also ignore the vacuum term, concentrating 
only on the matter contribution. Noting the symmetry in the matrix element under interchange of k and w, we may 
change variables w — > in the part of the integrand proportional to the distribution function n{w) to get 



8eV 
(2^)3 



dkh{k)Dq{k) 



{2TTf2q2k2w 



6{E 



2„2 



32e^(7 



E-2k E-2w 



){2n(fc)} 

E{E - 2q) 



E-2w E-2k [E-2w][E-2k] 



.17) 
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The argument of the delta function is the equation g(q) = k + q + w{q) — E = 0. The solution of this equation is 
at g = qs{k, E): 



The delta function can be written as 



1 E{E-2k) 
2E-k{l-y)- 



(8.18) 



Sigiq)) 



\9'{Qs)\ 



Substituting this back into the equation for Dq, we can do the dq integration with the above mentioned delta 
function. We can then perform the remaining angular integration by removing the small part of phase space e to 
shield the coUinear singularities. Now expanding up to linear order in k as k << E, we get 



2„2 



8e'g 



dkfi{k)Dq{k) 



8e^ 



dkn{k) 



2k 



2/c - £; ) log ( - 



Thus we get 



Dq{w) ^2w + (^~2w ~ E^ log (- 



(8.19) 



(8.20) 



We can now obtain Dg by concentrating on region d) of phase space and ignoring all combinations of distribution 
functions containing n{k) or h{w), this gives us, 



Dg{w) = -2w + ( 2w - 2^; + — ) log I - 



.21) 



D. Pair annihilation F & Compton scattering C 



The procedure to obtain these is almost exactly identical to the two terms of the previous section. The total Compton 
scattering cont ribut ion can be obtained from Eq. (6.17) by doubling it as mentioned in the paragraph immediately 
following Eq. ( 3.17 ). We may, once again, from phase space considerations show that the dominant contribution to 
Compton scattering occurs from a region where fc ~ T << .^(k is the incoming quark or antiquark energy). The 
leading term of Compton scattering is, thus, proportional to the quark or antiquark distribution function. From 
similar considerations the leading term of pair annihilation can be demonstrated to be proportional to the outgoing 
gluon distribution function. Expanding them up to linear order in the quark or gluon energy we get: 



C{w) ^2w+ (^-2w + E^ log 0^ 



.22) 



and 



F{w) = -2w + {2w + 2E 




(8.23) 



IX. RESULTS 



In the previous seven sections we evaluated the two different self-energies of the photon at two loops; then evaluated 
the various cuts of the self-energies which constituted its imaginary part; we then recombined the various cuts and 
reinterpreted them as physical processes; finally we evaluated these terms in the limit of heavy photon emission. In 
the last section we have concentrated solely on the thermal or matter part of these expressions. The vacuum part is 
well known. All the expressions contain coUinear singularities(as e — > 0), which for the moment have been shielded 
by removing the small part of phase space(e) where these singularities occur. Some of the expressions also display 
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infrared singularities as it; — > 0. Hence, the final integrations over w are yet to be performed. In the following we will 

combine all these terms and perform this integration. 

We now re-substitute the terms F, C, D, V, and S back in Eq. (8.1) to get the coefficients of the bosonic and 
fermionic distribution functions as 



n{w)(^F + Dg + Vg + Sg'^ =n(u;)(-4w). 



(9.1) 



liw) (C + Dq + Vg + Sg^ = n{w){~~8w). 



Thus, we find that when all the cuts are summed, the coUinear and infrared singularities cancel, 
contradistinction with the results of Refs. (|,|^, where the infrared singularities cancel but the coUinear 
persist. With these, we get the full imaginary part of the self-energy as 



(9.2) 

This is in 
singularities 



2^2 



ImH^ 

^ 21oop, thermal 



4e/g 
" 3 ■ 

We may also derive the Born term and quote the two-loop vacuum contribution(from Jl7[) as 

-3p2 

ImH'^ + ImH'^ 

'^lloop '^2 loop, vacuum 




(9.3) 



(9.4) 



X. DISCUSSIONS AND CONCLUSIONS 



In this paper, we have calculated the imaginary part of the two-loop heavy boson retarded self-energy in the 
imaginary time formalism. We also elucidate the analytic structure of the self-energy by recombining and reinterpreting 
various cuts of the self-energy as physical processes which occur in the the medium. Cuts with loops have been 
interpreted as interference terms between 0{a) tree scattering amplitudes and the Born term with spectators. At 
each stage the results from the self-energy cuts was matched by re-deriving the amplitudes of the tree level diagrams. 
This constitutes an important check of each of the contributions from the self-energy. 

Each of the contributions contain infrared and coUinear singularities. In the interest of simplicity, we analyzed 
this singular behaviour in the region where the dilcpton mass is far greater than the temperature. This allows us 
to neglect a series of terms which appear sub-dominant. In each case we retained terms only up to order T'^/E'^. 
One might argue that this represents a considerable approximation of the result. However the resulting simplification 
allows us to analyze far simpler and analytically integrable expressions. We would point out that this was precisely 
the approximation used in |^,^ where a remnant coUinear divergence was deduced at 0{T^ /E'^). When all the 
contributions were summed, all infrared and coUinear divergences cancelled; leaving a finite result of 0{T'^ / E'^). In 
this sense our results differ slightly from those of Ref. |^] who find a remnant 0{T^/E^) result. The possible reasons for 
this discrepancy are many. For example, the authors of Ref. use a complicated finite temperature remormalization 
prescription, ours is the same prescription as at zero temperature; they apply finite self-energy corrections on the 
outer legs of all their processes and we do not. However, both calculations (as well as those of Ref. [||) yield the 
consistent result that all coUinear and infrared divergences cancel in the final rate expression. This is consistent with 
the KLN theorem ||l|,|l|l, even though a formal proof of the theorem at finite temperature is still elusive. We leave 
this and other extensions for future investigations. 
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XII. APPENDIX A 



Notation 

Our notation is categorized by the explicit presence of an apparent Minkowski time = —ir and a momentum 
q° = i2mrT or i{2n + 1)ttT for bosons or fermions respectively. Our metric is (1, —1, —1, —1). For the case of zero 
chemical potential our bosonic propagators have the same appearance as at zero temperature, i.e., 

The Feynman rules are also the same as at zero temperature, with the understanding that we replace the zeroth 
component of the momentum by i(2n + 1)ttT for a fermion and by an even frequency in the case of a boson. One 
may, in the case of zero chemical potential, relate this to the familiar case of reference by noting that 

where A e {^n , q) is the familiar Euclidean propagator presented in the literature ( [p^ , |l5| ) . One may immediately 
surmise the form of the non-covariant propagator A(|g|, the Fourier transform of which is the covariant propagator. 

JO 

dx"e^9°-°A(|q1,a;"). (12.3) 

"'0 

The full fermionic propagators are 

^(9,9°) = (7%) (-0 / "'da;Oe''?"-°A^(|q1,xO), (12.4) 

JO 

where, 

A^(k1,:r°) = -^^/,(iJje--"(^^). (12.5) 



XIII. APPENDIX B 

Discontinuity across a second order pole. 

Imagine we have a function of a complex variable F{z), and it is given to be in the form 

Fiz)^ fdxl^ + p^, (13.1) 

J Z — X (z — X) 

where x is a real variable, integrated on the real axis. Most of the discontinuities that we evaluate can be cast in 
this general form. This can be rewritten as 

s f , f(z,x) q(z,x) 
- F{z) = / dx^^^ - p ' 13.2 

J X — z [x ~ zy 

The functions f{z,x) and g{z,x) are analytic in x and hence admit a Taylor expansion. 

df Id"^ f 

f{z,x) = f{z,x^ z) + -^{z^x = z)[x ~ z\ + 2^(^'^ ^ ^)[^ ~ + ■■■ {l^.?>) 



Substituting Eq.(13.3) in Eq.(13.2) we get 
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F{z) 



f{z,x — z) df d'^f , NT 1 

dx \- —(z,x = z) + -r^lz, X = z)\x - z\ + ... 

X — z dx dx^ 

1 A. _ , 1 

2dx^ ^'"''"'^'{x-zf 



g{z, X = z) dg 1 
[x — zy dx X ~ z 



(13.4) 



Recalling that only the pure first order poles develop a discontinuity or imaginary part at the pole we get the 
imaginary part of Eq.(13.4) as 



f{z,x) - ^iz,x) 
dx 



(13.5) 
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